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Resume 

We consider singular perturbations of elliptic systems depending 
on a parameter e such that, for e = the boundary conditions are not 
adapted to the equation (they do not satisfy the Shapiro - Lopatinskii 
condition) . The limit holds only in very abstract spaces out of distribu- 
tion theory involving complexification and non-local phenomena. This 
system appears in the thin shell theory when the middle surface is el- 
liptic and the shell is fixed on a part of the boundary and free on the 
rest. We use a heuristic reasoning applying some simplifications which 
allow to reduce the original problem in a domain to another problem 
on its boundary. The novelty of this work is that we consider systems 
of partial differential equations while in our previous work we were 
dealing with single equations. 

1 Introduction 

This paper is devoted to a very singular kind of perturbation pro- 
blems arising in thin shell theory. Up to our knowledge, it is disjoint of 
relevant and well known contributions of V. Mazya on perturbation of 
domains and multistructures for elliptic problems including the Navier 
- Stokes system ([l^, [U], [13]), as the pathological feature of our pro- 
blem is concerned with ill-posedness of the limit problem, generating 
singularities out of the distribution space. So, it may be considered as a 
contribution to enlarge perturbation theory of Mazya. More precisely, 
the main purpose of this paper is to generaHze the previous work done 
on equations, see [7], [H] to systems of partial differential equations. 
The motivation for studying that kind of problems comes from the 
shell theory. It appears that when the middle surface is elliptic (both 
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principal curvatures have same sign) and is fixed on a part Fq of the 
boundary and free on the rest Fi, the "limit problem" as the thickness 
£ tends to zero is elliptic, with boundary conditions satisfying Shapiro 
- Lopatinskii (SL hereafter) on Fq but not satisfying it on Fi. In other 
words, the "limit problem" for e = is highly ill-posed. This patholo- 
gical behavior arises only as £ = 0. In fact, for £ > the problem is 
"classical". 

In such kind of situations, the limit problem has no solution within 
classical theory of partial differential equations, which uses distribu- 
tion theory. It is sometimes possible to prove the convergence of the 
solutions towards some Hmit m", but this "Hmit solution" and the 
topology of the convergence are concerned with abstract spaces not 
included in the distribution space. 

The variational problem we are interested in is : 

J Find u" eV such that, <E V . . 

\ a{u'',v)+e%iu',v) = {f,v), ^ ' 

or, equivalently, the minimization in V of the functional 

a(u, u) + e'^b{u, u) — 2{f, u), 

where f G V is given and the brackets denote the duality between V' 
and V. 

This is the Koiter model of shells, £ denoting the relative thickness. 
The corresponding energy space F is a classical Sobolev space. 

The limit boundary partial differential system associated with l|l.ip 
when £ = is elliptic and ill-posed. 

Let us consider formally the variationnal problem of the membrane 
problem (i.e. £ = 0) : 

J Find u &Va such that, Vw G Va ,^ „x 

{ a{u,v)^{f,v), ^^-^^ 

where Va is the abstract completion of the "Koiter space" V with the 
norm |jw||a — a{v, vy^^, it is to be noted that the elements of T4 are not 
necessarly distributions. The term "sensitive" originates from the fact 
that this latter problem is unstable. Very small and smooth variations 
of / (even in ^{n)) induce modifications of the solution which are large 
and singular (out of the distribution space). 

The plan of the article is as follows. After recalling the Koiter 
shell model (Section [2)), we recall the definitions of ellipticity and 
the Schapiro-Lopatinskii condition for systems eUiptic in the DougHs- 
Nirenberg sense (Section [3]) . In Section [H we study four systems of 
partial differential equations which are involved in our study of shell 
theory. These systems are the rigidity system, the membrane tension 
system, the membrane system and the Koiter shell system. 

In section [5l we study a sensitive perturbation problem arising in 
Koiter linear shell theory and we briefly recall some abstract conver- 
gence results. In Section [U we report the heuristic procedure of [7J. 
In this latter article, we addressed a model problem including a va- 
riational structure, somewhat analogous to the shell problem studied 
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here, but simpler, as concerning an equation instead of a system. It is 
shown that the hmit problem involves in particular an elhptic Cauchy 
problem. This problem was handled in both a rigorous (very abstract) 
framework and using a heuristic procedure for exhibiting the structure 
of the solutions with very small e. The reasons why the solution goes 
out of the distibution space as e goes to are then evident. The heu- 
ristic procedure is very much analogous to the method of construction 
of a parametrix in elliptic problems [2^ , [8] : 

-Only principal (with higher differentiation order) terms are taken 
into account. 

-Locally, the coefficients are considered to be constant, their values 
being frozen at the corresponding points. 

-After Fourier transform (x ^), terms with small ^ are neglected 
with respect to those with larger ^ (which amounts to taking into 
account singular parts of the solutions while neglecting smoother ones) . 
We note that this approximation, aside with the two previous ones, lead 
to some kind of "local Fourier transform" which we shall use freely in 
the sequel. 

Another important feature of the heuristics is a previous drastic 
restriction of the space where the variational problem is handled. In 
order to search for the minimum of energy, we only take into account 
functions such that the energy of the hmit problem is very smah. This 
is done using a boundary layer method within the previous approxima- 
tions, i.e. for large |^|. This leads to an approximate simpler formulation 
of the problem for small e, where it is apparent that the limit problem 
involves a smoothing operator and cannot have a solution within dis- 
tribution theory. 

Notations are standard. We denote 

dk^^, fc = l,2, (1.3) 

and 

d 

Dk^-ijT—, k = 1,2 and D" = D"^D^\ a ^ {ai,a2) eZl. (1.4) 

Moreover, the definition of the Sobolev space H^{T), s e M, where T 
is a one dimensional compact manifold is classical using a partition of 
unity and local mappings. 

The inner product and the duality products associated with a space 
V and its dual V will be denoted by (., .) and (., .) respectively. 

The usual convention of summation of repeated indices is used. 
Greek and latin indices will belong to the sets {1, 2} and {1, 2, 3} res- 
pectively. 

2 Generalities on the Koiter shell model 

Let ft he a bounded open set of with smooth boundary P. Let 8^ 
be the euclidean space referred to the othonormal frame (O, ei, 62, 63). 
We consider the shell theory in the framework of the Koiter theory and 
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more precisely the mathematical framework of this linear theory. The 
middle surface S of the shell is the image in 8^ of for the map 

The two tangent vectors of S at any point y are given by : 

Oq = da^p, a g {1,2}, 

where da denotes the differentiation with respect to y" , while the unit 
normal vector is : 

ai A 02 

as = T 7 u- 

\\ai A a2\\ 

For simplicity, we omitted y in the previous notation {aa{y)). 

The middle surface S is assumed to be smooth (C°°) and we may 
consider in a neighbourhood of it a system of "normal coordinates" 
y^,y'^,y^, when y^ is the normal distance to S. More precisely we consi- 
der a shell of constant thickness e, i.e. it is the set 

C^{Me£\ M = ipiy\y^)+y^a3, {y\y^) e il, -^e < y"" < ^e}. 

Under these conditions, let u = u{y^,y'^) be the displacement vector 
of the middle surface of the shell. In the linear theory of shells, which 
is our framework here, the displacement vector is assumed to describe 
the first order term of the mathematical expression as the thickness e 
is small, see [4} 118). 

Remark 1. In the sequel smooth should be understood in the sense of 
C°°. 

Remark 2. We consider here the case where the surface is defined by 
only one chart but this could be easily generalized to the case of several 
charts (atlas). 

More precisely, since we consider the case where u is supposed to 
be small, the Koiter theory is described in terms of the deformation 
tensor (or strain tensor) jaf3 of the middle surface : 

7a/3 — 2 ~ ""/3) 

and the change of curvature tensor pap '■ 

Pap = bap — bap- 

In the previous definitions, the expressions Uap (resp. dap) denote the 
coefficients of the first fundamental form of the middle surface before 
(resp. after) deformation : 

aap = aa- ap = daP ■ dpp, 

and bap (resp. bap) the coefficients of the second fundamental form 
accounting for the curvatures before (resp. after) deformation : 

bap = -Oq • dpa^ = as • dpUa = as • daap = bpa, 



4 



due to the fact that • as = 0. 
The dual basis a* is defined by 

ai-a^ =Sl, 

where 6 denotes the Kronccker symbol. The contravariant components 
a*-^ of the metric tensor are : 

a'^ = a' ■ a', 

and Uij are used to write covariant components of vectors and tensors 
in the usual way. Finally, the tensors 7 and p take the form : 

7/3a(u) = lafliu) = ^{Ua\fj + U^|a) - &a/3M3, (2.1) 
Paf3{u) = U3\oci3 + b0\„UX + b^Ux\a + b^Ux\/3-b^bxf3U3, (2.2) 

where S^as = b^^a^, b^ — (i^^h(xa^ .la denotes the covdricLTit differen- 
tiation which is defined by 



and 



U3\f3 = df)U3, 



bl\p = d<^b>^p + ^Lb''p-T^0y.=b^^^„ 

U3\afS = dapU3 - TKdxU3, 



(2.3) 
(2.4) 



where F^^ are the Christoffel symbols of the surface 
= F«^ = a« • dua^ = a« • d^af}. 

Let us now define the energy of the shell in the Koiter framework. 

It consists of two bilinear forms a and h : a corresponds to a membrane 
strain energy and 6 is a bending energy (which acts as a perturbation 
term). More precisely, a is defined by 

a{u,v)= [ A"''^''7A^(«)7aMt') ds, (2.5) 



where A"'^^^ are the membrane rigidity coefficients which we assume 
to be smooth on il. Moreover, we assume that some symmetry holds 

Defining the membrane stress tensors by 

r"^'(u) = A"^^^7a4u), (2.7) 
using the symmetry of 7, we immediately see that 

= T'^"{u), (2.8) 

and 

a{u,v) = [ T"^{u)jc.0{v) ds= [ 7a/3(w)T"'3(t7) ds. (2.9) 



Furthermore, we assume that a coercivity condition holds uniformly 
on the surface : 

^"''^''^a/36/.>C||$f , C>0. (2.10) 
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Remark 3. It is to be noticed that there are two different symmetries 
on A : the first one A"^^^ — A^^^"-'^ is necessary to exchange u and v 
in (122]) while the second A^^"^'^ = A^^"^ is used to obtain but is 
not necessary in order to obtain l|2.9p since we could use the symmetry 
of 7. 

Analogously, we define the bilinear form b which corresponds to the 
bending energy of the shell and which will act as a perturbation term : 

b{u,v)= f B"P^''p^^{u)po.p(v) ds, (2.11) 



where B"^^^ are the bending rigidity coefficients which we assume to 
be smooth on Q and to have the same properties l|2.6p and l|2.10p as 
A, namely 

j^apXji _ j^\fiaf3 _ j^tiXafS ^2 12) 

and 

B'^^^^^c.pix^ > CUW^ (2.13) 

uniformly on the surface. 

Similarly to a we can write 

b{u,v)^ [ M''^{u)p^p{v)ds, (2.14) 
Js 

where the bending stress tensors are 

M"^(u) = S"^^^pam(w)- (2.15) 

In this work, we will restrict ourselves to the case of elliptic surface, 
i.e. we will always assume that the coefficients bai3 are such that 

^11^22 ~ &12 > uniformly on S and 611 > 0. (2-16) 

Let us finish this introduction by topoligical considerations, the 
boundary dO, = Fq U Fi is assumed to be smooth (i.e. of class C°°) 
in the variable y = {y^,y'^), where Fq and Fi are disjoint; they are 
one-dimensional compact smooth manifolds without boundary, then 
diffeomorphic to the unit circle. 

We consider the following variational problem (which has possibly 
only a formal sense) 



Find u"" eV such that, Vv e V 

a{u'',v) +e^b{u^,v) = {f,v), 



2Un,e n,\ - / f „\ (2-17) 



with a and b defined by (|2.5p and l|2.14p where the space V is the 
"energy space" with the essential boundary conditions on Fq 

V = {v; Va £ H^{Q), V3 e H^{Q); v^^o = in the sense of trace}. 

(2.18) 
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Remark 4. The essential boundary conditions on Fq l|2.18p corresponds 
to the case of the fixed boundary of the shell. Other boundary condi- 
tions could have been considered such as : 

V = {v; Va e W3 £ H^{Q); viro = 0, d^v^^^o = in the sense of trace}, 

(2.19) 

where v is the normal to Fq (i.e. the normal to the boundary which 
lies in the tangent plane), which corresponds to the clamped case. 
The following Lemma was obtained by Bernardou and Ciarlet see 

Lemma 2.1. The bilinear form a + b is coercive on V. 

We shall denote by V the dual space of V. Here dual is obviously 
understood in the abstract sense of the space of continuous linear func- 
tional on y . In order to make explicit computations in terms of equa- 
tion and boundary conditions, we shall often take / as a "function" 
defined on fi, in the space 

{/ e H-^n;R) X H-i(f7;K) x H-^{n;R); (2.20) 
/ "smooth" in a neighbourhood of Fi} C V , 

where "smooth" means allowing classical integration by parts. Ob- 
viously other choices for / are possible. 

Moreover, we immediately obtain the following result. 

Proposition 2.2. For s > and for f in V , the variational problem 
\2.1'7\) is of Lax-Milgram type and it is a self-adjoint problem which 
has a coerciveness constant larger than ce^ , with c > 0. 

Remark 5. It is to be noticed that the coerciveness of the previous 
problem disapears when e = 0. 

3 The ellipticity of systems and the Shapiro- 
Lopatinskii condition 

In this section, we recall some classical results on the linear boun- 
dary value problems for elliptic systems in the sense of Douglis and 
Nirenberg We begin with the definition of ellipticity for systems, 
then we recall the Shapiro-Lopatinskii condition. This latter condition 
states which boundary conditions are well suited in order to have well 
posed problems for elliptic systems. We then recall in what sense an 
elliptic system with Shapiro-Lopatinskii condition is "well-behaved". 

For brevity, from now on we will denote SL the Shapiro-Lopatinskii 
condition. 

3.1 Elliptic systems in the sense of Douglis and Ni- 
renberg ^ 

In this work, we shall deal with systems of ^ (Z = 3 or Z = 6) 
equations with 3 unknowns (noted here mi, U2, U3) defined on an open 
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set O C with smooth boundary, which has the form : 

hjUj ^ fk, k ^1,...J, (3.1) 

or equivently Lu — f. The coefficients lkj(x,D) with D = {Di,D2) 
and Di = —i^-, I G {1,2}, are Unear differential operators with real 
smooth coefficients. In our systems p.ip . the highest order of differen- 
tiation is different for the three unknowns and depends on the equation. 
A way to take into account such differences between the various equa- 
tions and unknowns is to define integer indices (si, S2, S3) attached to 
the equations and integer indices (^1,^2,^3) attached to the unknowns 
(see Doughs and Nirenberg [6J) so that the "higher order terms" (which 
will be called "principal terms") are in equation j the terms where each 
unknown "k" appears by its derivative of order +tj. More precisely, 
the integers {sk,tj) are such that 

J if Sfc -f tj > 0, the order of Ikj is less or equal to Sk + tj, 
\ if Sfc + tj < 0, Ikj is equal to zero. 

The principal part 1'f.j of Ikj is obtained by keeping the terms of order 
Sk + tj if Sk + tj > and by taking 1'f.j = if Sfc -I- tj < 0. The matrix 
L'{x,S), ^ — (^1, ^2) G R^, obtained by substituting for in ZJ,^ , is 
called the principal symbol of the system. Since I'f^j are homogeneous of 
order Sfc-|-tj with respect to ^a, the determinant of the matrix L'[x,£), 
denoted D(x,^), is homogeneous of degree SfeSfc + Tijtj. 

Definition 3.1. The system l|3.ip is elliptic in the sense of Douglis 
and Nirenberg at the point x G fl ii and only if 

B{x,O^0, VCeK'\{0}. (3.2) 

Remark 6. Since the coefficients are assumed to be real, the function 
D(a;,^) for an elliptic system is even in ^ of order 2m with 

SfcSfc -f Sjtj — 2m. 

Remark 7. The definition of the indices Sj and tk for a system is 
slightly ambiguous. Indeed the result is exactly the same after adding 
an integer n to the indices Sj and substracting n from the tk- 
Remark 8. Let xo e be such that the system p.ip is not elliptic, 
then there exists a ^ G \ {0} such that D(a;o,<^) = 0. In such a 
case the system L'{xq,D)u = 0, with frozen coefficients at Xq admits 
a solution of the form u{x) = we*^^, with u e M'^ \ {0}. 
Remark 9. Moreover, throughout this paper, ellipticity will be unders- 
tood in the sequel as uniform, i.e. there exists a positive constant A 
such that 

^"'SalCaP < |deti'(x,C)| < AE.ICcP, 
for &\\x&n and ^ = (Ci, 6) e 
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3.2 Shapiro-Lopatinskii conditions for elliptic sys- 
tems in the sense of Douglis and Nirenberg [6j 

From now on, for simplicity, we will say that a system is elliptic 
when it is elliptic in the sense of Douglis and Nirenberg [6J. 

Let Ikj (L) be an elliptic system of order 2m with principal part I'f.^ 
{L') and let m boundary conditions be given by : 



where bkj {x, D) are differential operators with smooth coefficients. 
Let us define the integers (indices of the boundary conditions, 
k = 1, .., m) such that 



The principal part h'^.^ is h^j if +tj>Q and zero otherwise. 

Assume that the smooth real coefficients are defined in fi. 

Let xo £ r, we assume that L' is elliptic at xq. Usually, see and 
[H] for instance, the SL condition at a;o is defined via a local diffeom- 
rophism sending a neighbourhood of xq in fl into a neighbourhood of 
the origin in a half-plane. For ulterior computations, it is worth-while 
to take a special diffeomorphism which amounts to taking locally car- 
tesian coordinates xi,X2, respectively, tangent and (inwards) normal 
to the boundary at xq. We then consider only the principal parts of 
the equations and of the boundary conditions frozen at a;o. Next, we 
consider the corresponding boundary value problem obtained by formal 
tangential Fourier transform (i.e. Di ^ ^i, with S M and u ^ u) 
which amounts to the following algebraic conditions : 



j, k e {1, m}, see [7J Sec. 3.2 for details, if necessary. 

The problem (|3.3p involves a system of ordinary differential equa- 
tions with constant coefficients of the variable X2 G M"'" and m boun- 
dary conditions at X2 = 0, whose solutions are classically a linear 
combination of terms of the form : 



P(x2)e*^^^^, where P is a polynomial, in the case of Jordan block. 



Recalling that the system L is elliptic, it follows that the imaginary 
part of ^2 does not vanish. Furthermore, there are m solutions ^2 of 
D{xq,S,i,£^2) = with positive imaginary part that we denote (and 
m with negative imaginary part denoted ^^). 

We then try to solve l|3.3p using only linear combinations of the 
m solutions of the form p.4p for the m roots (i.e. exponentially 
decreasing towards the domain). 



hkjUj = gk, k e {1, ..,m}. 



if rk + tj > 0, the order of bkj is less or equal to -I- tj 
if Tk +tj < 0, bkj is equal to zero. 




(3.3) 




(3.4) 
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Definition 3.2. The SL condition is satisfied aX xq G T if one of the 
following equivalent conditions holds : 

1. The solution of the previous problem is defined uniquely. 

2. Zero is the only solution of the homogeneous (i.e. with gj = 0) 
previous problem. 

Remark 10. The two conditions (which are equivalent) of the previous 
definition are clearly equivalent to the non annulation of the determi- 
nant of the corresponding algebraic "system". 

Remark 11. The reason for defining the SL condition amounts to the 
possibilty of solving the problem in a half plane via tangential Fourier 
transform. The reason for not considering the roots is that, for 
X2 > 0, they should give exponentially growing Fourier transforms in 
xi ^1, which are not allowed in distribution theory (note that 
and ^2 are proportional as D{^) is homogeneous). 

The verification of the SL condition is often tricky. In some situa- 
tions, we can use equivalent expressions which are simpler to treat. 
More precisely, define the function u by u{xi,X2) = u(^i, a;2)e'I^^I^S 
with u{^i,X2) = ve^^^ (or expressed as exponential polynomial in the 
case of Jordan block) , it is an exponentially decreasing function in the 
direction inwards the domain (when X2 —>■ +oo), it is also a periodic 
function in the tangential direction Xi and it satisfies 



j,k G {l,...,m}. The following proposition is very useful in the case 
where ellipticity is linked with positive energy integrals obtained by 
integrating by parts. For instance, we have : 

Proposition 3.3. Consider the homogeneous problem associated with 
f3. 5\) (i.e. taking gj — 0) for xq G T. If any solution u, which is 
periodic in the tangential direction xi and exponentially decreasing in 
the direction X2 inwards the domain, is zero, then the SL condition is 
satisfied. 

Remark 12. In order to have well-posed problems for elliptic systems, 
boundary conditions satisfying the SL condition should be prescribed 
at any points of the boundary. Their number is half the total order of 
the system. 

Remark 13. The specific boundary conditions may differ from a point 
to another on the boundary. In particular, each connected component 
of the boundary may have its own set of boundary conditions. Other- 
wise, local changes of boundary conditions (as well as non-smoothness 
of the boundary) induces local singularities. A changer 

3.3 Some results for "well posed" elliptic systems 

Let us now consider a boundary value problem formed by an el- 
liptic system with boundary conditions satisfying the SL condition. In 
what sense is this problem "well-behaved" ? The obvious example of 




(3.5) 
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an eigenvalue problem, even for an equation shows that uniqueness 
is only ensured up to the kernel formed by the eigenvectors associa- 
ted with the zero eigenvalue, whereas existence involves compatibility 
conditions (orthogonality to the kernel of the adjoint problem). The 
general results are those of Agmon, Doughs and Nirenberg [l]. 
First, let us recall the definition of a Fredholm operator. 

Definition 3.4. Let E and F be two Hilbert spaces and A an operator 
(closed with dense domain in E) from E into F. We say that A is a 
Fredholm operator if and only if the following three conditions hold : 

1. Ker(A) is of finite dimension, 

2. R{A) is closed, 

3. R{A) is of finite codimension. 

The operator A is also said to be an index operator, the index is defined 
as dim Ker(A) — codim R(A). 

Let us consider an eUiptic system of order 2m whose coefficients 
are smooth : 



whose indices associated with unknowns, equations and boundary condi- 
tions are respectively tj,Sj,rj. Let p be a "big enough" real number, 
called regularity index. Consider operator p.6p as a linear operator 
from the space E to the space F defined by : 

E = ^^^HP+^^{n), F = ^^j^Hp-'^in) X irjij^Hp-'''-^dn). (3.7) 

The real p is chosen in order to give a sense to the traces which are 
involved, i.e. it is such that p — rj — 1/2 > for j G {!,..., m}. 

The following result is the main result of the theory of Agmon, 
Doughs and Nirenberg : 

Theorem 3.5 (Agmon, Doughs and Nirenberg [1]). Let fl be a bounded 
open set with smooth boundary T. Let us consider an elliptic system 
with boundary conditions satisfying the SL condition everywhere on T. 
Assume that the coefficients of the system are smooth and that u, f and 
g satisfy i3. 6\) . Then the following estimate holds true : 



where C does not depend on u,f,g. Moreover, the operator defined by 
lis. 6]} from the space E to the space F, given by iS. ?)) ■ is a Fredholm 
operator, for all value of p such that p — rj — 1/2 > for j G {1, rn}. 
Furthermore, the dimension of the kernel and the dimension of the 
subspace orthogonal to the range do not depend on p. The kernel is 
composed of smooth functions. 

Remark 14. The previous theorem means that in general existence 
and uniqueness of the solution only hold up to a finite number of 
compatibility conditions for / and g and existence of the solution holds 
up to a finite dimension kernel. More precise properties need specific 
properties of the system. 



{ 



IkjUj = /fe, j, k e {1, ...,1} in n 
bhjUj ^ gh, he {1, m} on dil, 



(3.6) 



\\u\\e < C{\\{f,g)\\F + ||w||(i2(o))i), 



(3.8) 
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Remark 15. For all values of p, the kernel formed by the eigenvectors 

corresponding to the eigenvalue is of finite dimension and is composed 

of smooth functions, independent of p (in C°°{fl)). 

Remark 16. Denote A the operator defined by l|3.6p in the spaces E 

and F. Let us consider the case where dim Ker(A) > and define the 

inverse i? of A as a closed operator from R{A) to £'/Ker(A), we have 

that 

ll"IU/Ker(A) <C||(/,5)||f, (3.9) 

where u is an element of the equivalence class of u. 

The element u can also be viewed as an element of the orthogonal 
of Ker(A) in E, which is identified with E/KeT{A). In such a case, 
there exists a unique {u,u) G E/KeT{A) x Ker(A) such that 

u = it + u. 

Since Ker(A) is of finite dimension, all the norms are equivalent and 
we can choose for ii a norm in a space H~'^ with v very big. Therefore, 
inequality p.Sp can be rewritten as 

\\u\\e < C{\\u\\E/KeT(A) + \\u\\H—) < C:{\\u\\ e /KeT{A) + \\u\\ H—) , (3.10) 

for v big enough such that E C H^". Recalling l|3.9p . we then deduce 
that 

\\u\\E<C{\\if,g)\\F + \\u\\H-^)- (3.11) 

Moreover, the norm in H^'^ may be replaced by a seminorm, provided 
it is a norm on Ker{A). 

Remark 17. In the case where dim Ker(yl) — 0, the inverse B of the 
operator A is well defined on R{A). It is a closed operator, hence it is 
bounded and the following estimate holds : 

\\u\\E<C\\{f,g)\\F. (3.12) 

4 Study of four systems involved in shell 
theory 

In this section, we study four systems, denoted by rigidity system, 
membrane tension system, membrane system and Koiter shell system, 
which will appear in the sequel. We prove that these four systems sa- 
tisfy the ellipticity condition and we study some boundary conditions. 
It is to be noticed that the boundary conditions may be different on 
Fq and Fi which are supposed to be disjoints. 

Let us recall the situation : is a connected bounded open set of 
R2 with C°° boundary F = FgUFi and FonFi = 0. The middle surface 

5 of the shell is the image in £^ of fl for the map 

We assume that the ellipticity assumption of the surface holds : 

~ ^12 > uniformly on il. 
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4.1 The rigidity system 

Let us begin with the rigidity system defined by ^ap(u) : 



Clearly Ua and U3 play very difi^erent roles as Uq, appears with deri- 
vatives whereas U3 only appears without. Therefore take (1,1,0) as 
the indices of the unknowns (mi,M2,U3) and (0,0,0) as equation in- 
dices in the order (711,722,712)- The principal system is obtained by 
substituting for F"^ but keeping 5a^. 

Lemma 4.1. Do to the ellipticity assumption of the surface 112. 16\} . 
the rigidity system 7 is elliptic of total order 2 on fl. 

Demonstration. Substitute —i£,a for da in the principal system, we 
obtain a system whose determinant is D{x, £,) = 2&i2^i$2— ^22^1 — ^ii^l , 
hence due to the ellipticity hypothesis l|2.16p . for all a; g il, we have 



4.1.1 Cauchy boundary conditions 

It is classical that the Cauchy problem associated with elliptic sys- 
tem is not well posed in the sense that it does not enjoy existence, uni- 
queness and stability of solutions. Nevertheless, the Cauchy problem 
associated with the rigidity system will be involved in the sequel and 
we study it now. In particular, we shall need the following uniqueness 
theorem for solutions u G x x L^. 

Lemma 4.2. Under the ellipticity assumption of the surface !12.16\} . 
the system 7q/3(u) — on fl with the boundary conditions ui — U2 — 
on apart of the boundary (of positive measure) admits a unique solution 
which is u — 0. 

Demonstration. Let us assume that v £ {il) x (fi) x H'^ (fi) is such 
that 7q/3(i') = and wi = f2 = on a part of the boundary. Thanks 
to the ellipticity hypothesis (|2.16p . we know that 611 7^ on fl. We 
can eliminate V3 from the first and third equations (711(1;) = and 
722(1') = 0) of the system 7. This yields the system of two equations 
for two unknowns {vi,V2) '■ 




(4.1) 



D{x,O>0- 



□ 



( 



= d2V2 - r^2V» - ^77(51^^1 - ^ll^c) 

= l{d2Vi + diV2) - T^2Vo. - j^idivi - T'^^va). 



(4.2) 



The eliminated unknown being then given by : 



V3 = 7—{dlVl - r"iWa). 
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The problem then reduces to the uniqueness in the class H^{Q) of 
{vi,V2) satisfying 

diVi — biiv^j ~ 

d2V2 - = (4.3) 

^{d2Vi + diV2) - bi2V3 = 0, 

with vi = W2 = on a part of the boundary. This problem is more or 
less classical. Under analyticity hypotheses about the coefficients and 
the boundary, the uniqueness follows from Holmgren local uniqueness 
theorem and analytic continuation (as mi,U2 are in this case analytic 
inside ft). Under the C°° hypotheses adopted here, uniqueness follows 
from theory of pseudo-analytic functions. There are two nearly equi- 
valent theories of such functions attached to the names of L. Bers (see 
for instance supplement of chapter IV of [5] , written by Bers himself) 
and I.N. Vekua see pT] . 

Let (^1,^2) be a solution of (|4.3p vanishing on a part T of the 
boundary. Let {vi,V2) be an extension of {vi,V2) with values zero to 
an extended domain across F. Classically {vi,V2) satisfies the same 
system (|4.3p on the extended domain and, according to interior regu- 
larity theory for elliptic systems, is of class C°° inside it. The function 
w = vi + iv2 is pseudo-analytic, of class C°° and vanishes on the ou- 
ter region of the extended domain. We then use either theorem 3.5 of 
|21) . p. 146, which gives directly the uniqueness or the representation 
theorem of [5] p. 379. In this case, w(z) admits the expression (here 
z = Xi + 1x2) ■ 

w{z) = e^(^)/(z), 

where f{z) is analytic and S{z) is continuous. As e"^^^^ vanishes now- 
here, the uniqueness follows. □ 

Remark 18. Strictly speaking, the evoked theorems of pseudo-analytic 
functions apply to systems with principal part of the canonical form 

I f>'"'~A""'^'" (4-4) 

[ d2Vi + diV2 = 

so that the classical reduction to this form (see for instance [5] p. 
169-170) should be previously considered. But obviously, this does not 
modify the C°° regularity inside the domain. 

Let us make several comments about this uniqueness result. 
Remark 19. This result, known as the infinitesimal rigidity of the sur- 
face, does not depend on the curvilinear coordinates. 
Remark 20. The key ingredients of the previous uniqueness result are a 
uniqueness theorem for the Cauchy problem for elliptic systems of two 
equations of order 1. It is not based upon a coercivity assumption for 
an elliptic system. But we know that the Cauchy problem for elliptic 
systems is precarious in the sense that it does not enjoy existence, 
uniqueness and stability of solutions. This means that such a system 
could lead to instability in the sense that there could exist wi,W2,i'3 
very "big" in usual spaces such that 7q/3(w) are very "small". 
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4.1.2 Boundary value problems for the rigidity system 

From now on, we will consider the frame (0,ai,a2,a3) to be or- 
thonormal on the boundary and such that Ut = (ui,0,0) and u„ = 
(0,?X2,0), where Ut denotes the component of u in the tangential di- 
rection to the boundary and m„ is the component of u in the normal 
direction to the boundary and in the tangent plane. This point which 
is not absolutely necessary, implies a special local parametrization. 

Lemma 4.3. The boundary condition ui — g satisfies the SL condition 
for the system 7. 

Demonstration. We take as index of the boundary condition r = — 1. 
Let a;o belong to F. As explained in Section 13. 2[ using a partition 
of unity, local mappings, with axes yi tangential and j/2 inwards F, 
dropping lower order differential terms, we obtain a new system : 

{diUi - 611U3 = 
d2U2 - 622M3 = (4.5) 
i(92Ui + diU2) - bi2U3 = 0. 

We look for solutions which are exponentially decreasing when 1/2 
+00 of the form : 

u{yi,y2) = Ue'^y-+'i'y\ ^ieR\{0}, 

(u,\ 

with = U2 G C'', Im(C) > 0. Substituting this solution into 

l|4.5p and using the boundary condition we have Ui — 0. Consequently, 
ui = everywhere and l|4.5p gives also U2 = U3 = 0. U2 = ~ 0. □ 

Remark 21. Similarly to the proof of the previous result, we can prove 
that the following boundary conditions satisfy the SL condition : 

1. U2 = g (take r = —1). 

2. U3 = g (take r ~ 0). 

Remark 22. Since Fq and Fi are disjoints and thanks to the previous 
statements, the boundary value problem 

7q/3(u) = on n, 

U2 = on Fo, (4.6) 
U3 = {t on Fi. 

is "well posed" in the Agmon, Douglis and Nirenberg sense. Recalling 
Theorem 13. 51 and Remark [HI together with standard regularity theory 
for elliptic systems, it follows that u is of class C°° on O U Fq for 
any u (either smooth or not). Consequently, up to a kernel of finite 
dimension composed of smooth functions belonging to C°°{fl)^ (and 
eventually up to a compatibility condition (to belong to the range of 
the operator which is a closed subspace of finite codimension), the 
space {v, 7q/3(w) = on fi, Vn — on Fq} is isomorphic with the 
space C°°(Fi). The previous statements can be rephrased as follows : 
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up to a finite dimensional space composed of smooth functions, the 
space {v, 7q/3(w) = on il, w„ = on Fq} is isomorphic to the space 
of traces on Fi : 

{veC^iT^)}, (4.7) 

the isomorphism is obtained by solving l|4.5p . 

In the sequel, we shall consider indifferently the functions v (defined 
up to an additive element of the kernel) or their traces v on Fi. 

4.2 The system of membrane tensions 

Consider the membrane tensions system T of three equations with 
the three unknowns (T^^ ^T^"^ ,T'^^) : 

^^ll rp21 _ J-1 

"-'ll ^ ^|2 — J 

-T^^-T^^^=f (4.8) 

It is apparent that the three unknowns play analagous roles. Concer- 
ning the equations, it is clear that the first and the second are similar 
but different from the third. Therefore, we consider (1, 1,0) as indices 
of equations and (0,0,0) as indices of unknowns. The principal sys- 
tem Tp is obtained by replacing the covaraint differentiation |^ by the 
usual differentiation da (i.e. replacing F^^ by zero). Proceeding as in 
the proof of Lemma |4!T| we obtain the following result. 

Lemma 4.4. Under the elUpticity assumption of the surface f2.16]} . 
the system T is elliptic of total order two. 

Remark 23. It is worthwhile to study the Cauchy problem for the 
membrane tension system l|4.8p . This is done exactly as in Section 
4.1.2 for the rigidity system. We eliminate one of the unknowns, T^^ 
for instance and (|4.8p reduces to an elliptic system of two first order 
equations in T^^ and T^^. The Cauchy conditions are T^^ = T^^ = on 
a part of the boundary. According to our special frame, this amounts to 
T'^^Uf^ = 0. This Cauchy problem enjoys uniqueness but not existence 
and stability in usual spaces. 

Remark 24. The system of membrane tensions T (|4.8p and the system 
of rigidity 7 l|4.ip are adjoint to each other. This is easily checked by 
covariant integration by parts on S. Indeed, neglecting boundary terms 
(we are only interested in the equations) and using (|2.ip together with 
the symmetry of the T"^), we have : 

/ T^^Ja^iu) ds 

Js 




[ T{T)u ds 
Js 



16 



4.3 The membrane system 

We denote by membrane system the system of three equations with 
three unknowns u = (wi, U2, U3) obtained from (|4.8p when the tensions 
are written in terms of u, i.e. 

-T|f(u)-T2i(«)=/2 (4.9) 



with 
and 



r"^(u) = A"^^^7A^(«), (4.10) 



T|f (^) = d.T-f'iu) + rf„r""(^) + r^T^^iu). (4.ii) 



In order to prove the eUipticity of the membrane system, we replace it 
by another, equivalent one. Indeed, we shall take as unknowns Ui, U2, U3 
and the supplementary auxiliary unknowns T^^, T^^, T^^. Inverting the 
matrix A"'^^^^ in (|4.10p and recalling the definition of 7, we obtain the 
following equivalent system : 

^^ll 7^21 _ fl 
-'ll ^|2 —J 

-Tfi-Tf^=f (4.12) 
-6nr" - 2&i2ri2 - 622T22 = /3, 

W2I2 - b22U3 - C22a;3T"^ = (4.13) 
\{ui\2 + W2|l) - &I2W3 - C12Q/3T"'' = 0, 

where Capx^i are the compliances (inverse matrix oi A°'l^^^^). The sys- 
tem (|4.12p and l|4.13p is a system of six equations with the six unk- 
nowns (T^^, r^^, T^^, Ml, U2, U3) (written in this order). We recognize 
the membrane tension system in l|4.12p and the rigidity system in 
l|4.13p . Consider (1, 1, 0, 0, 0, 0) as indices of equations and (0, 0, 0, 1, 1, 0) 
as indices of unknowns. Then replacing the differentiation da by — 
and taking the determinant of the obtained system, we have a deter- 
minant of the form 



Dii 

£'21 -022 



= Dii LD2 



where the Dap are 3x3 matrices. Moreover, Dn and D22 are precisely 
those of the membrane tension system and the rigidity system respec- 
tively and ellipticity follows. The same result is obviously obtained 
without using the auxihary unknowns T"'^, in fact, we have. 

Lemma 4.5. Under the ellipticity assumption of the surface 112. 16\) . 
the membrane system with indices (of unknowns and of equations) 
(1,1,0) (1,1,0) is elliptic of total order four. 

Let us now state boundary value problems which will be considered 
later on. It is to be noticed that only two boundary conditions are 
considered on Tq. 
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Proposition 4.6. The boundary value problem 



< -buT^Hu) - 2foi2Ti2(u) - b22T^Hu) = f (4.14) 
ui — U2 = 0, on Tq 
T°'f'{u)na = on Ti, /? G {1, 2}. 

with unknown u satisfies the SL condition on Tq but it does not on Ti . 

Remark 25. The partial differential boundary value problem l|4.14p is 
formally associated with the variational problem (|2.17p when e = 0. 

Demonstration. Let us fix xo S F. According to the definition of the 
SL condition, we consider the homogeneous system with constant co- 
efficients in which we only kept the principal terms, i.e. taking F^^ = 
but bafs ^ and p = 0. 

After a change of coordinates with local mappings, still denoted by 
{xi,X2), we only have to consider solutions, which are exponentially 
decreasing in the direction inwards the domain (X2), of the correspon- 
ding boundary value problem obtained by formal tangential Fourier 
transform. Denoting by u{^i,X2) such a solution, by periodicity, we 
can restrict the domain to the strip B = (0, 27r/|^i|) x (0, +00) and we 
can consider the function 

v{xuX2)^e'^'^'u{^i,X2), (4.15) 

which is periodic in the tangential direction xi, decreasing as X2 +00 
and satisfies the homogeneous boundary condition associated with the 
principal part of (|4.14p . Recall that v satisfies the equation 

-(9iT"(w) -<92T21(w) = 

-a,T22(i;) - diT^^iv) = (4.16) 
-feiiT"(«) - 26l2Tl2(^;) - b22T^Hv) = 0. 

We multiply each line of l|4.16p by the conjugate Vi and we integrate by 
parts on the periodicity layer B. We see that on the infinite boundary 
the boundary integral is vanishing thanks to the decreasing condition as 
X2 +00. The boundary integral also vanishes on the lateral boundary 
(which is parallel to X2) of the strip thanks to the periodicity of v. 
Recalling the definition of T*^ , we obtain 

/ A"'5^^7A4z;)7a/3(^J) d:Ei dxa = 0, (4.17) 
Jb 

where obviously all the F^^ = 0. Consequently, recaUing the positivity 
property l|2.10p of A, this yields that 

/ Sa/3|7a/3(w)l' dxi da;2 = 0, (4.18) 
Jb 

and then 

Icfsiv) = on B. (4.19) 



18 



We have now to distinguish two cases. 

If xo e To, then reasonning as in Lemma 14.21 (or merely as in 
Lemma [4. 3p . we deduce that vi = V2 — V3 — 0, which means that the 
SL condition is satisfied on Fq. 

Let now a;o S Fi and 

7„^(u) = on B. (4.20) 

Remembering the definition l|4.15p of v, this yields that u is a solution 
of the following system of ODE of order 2 : 

i^iui - 611-U3 = 

92^2 - b22U3 = (4.21) 
i(92Wl + i6'"2) - bi2U3 = 0. 

Thanks to the fact that bu ^ and 622 7^ this can be rewritten as : 

biidlu2 - 21612^1 92^2 - &22C1W2 = 0. 

Recalling the ellipticity condition (|2.16p , we obtain after an easy com- 
putation that there exists a complex solution u, given by u = we^-^'^ , 
where w 7^ and A_ is the root with negative real part of 

6iiA2 - 2ibi2ii\ ~ 622^? = 0. 

This means that there exists non zero v which is exponentially decrea- 
sing in the direction inwards the domain 

v{£,i,X2) = W«i^ie^-^% 

with Re(A_) < such that 

lapiv) = on B, 

and hence 

T"^(w)n„ = on Fi. 
Therefore, the SL condition is not satisfied on Fi. □ 



4.4 The Koiter shell system 

The boundary value problem associated with the variational pro- 
blem l|2.17p with e > is classical and well-posed (see for instance [4], 
[TT]). It is eUiptic of total order 8, and the boundary conditions satisfy 
the SL condition. The system of equations is obtained by integration 
by parts, which yields : 




(4.22) 
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where the flection moments M°'P were defined in l|2.14p . (|2.15p The 
boundary conditions on Fq (supposed clamped) are : 

ui — U2 = M3 = dnU^ = on Fq (4.23) 

while the natural boundary conditions on Fi are in number of four, 
are not relevant (they are boundary terms obtained by integration by 
parts). We have : 

Proposition 4.7. The boundary value problem associated with the 
variational problem f2.17]} when e > considered as a system of three 
equations with the unknowns u is elliptic of total order 8 with indices 
(1, 1, 2) for the unknowns and the equations. 

5 A sensitive singular perturbation problem 
arising in the Koiter linear shell theory 

Very few is known concerning elliptic problems with boundary 
conditions not satisfying the SL condition and there is no general 
theory concerning them to our knowledge. Linear shell theory is one 
physical theory where they are naturally involved. 

5.1 Definition of the problem 

Let us first recall the variational problem l|2.17p we are interested 

in : 

f Find u" &V such that, 'iv , , 

\ a{u'',v)+e%{u'',v) ^ {f,v), > 

where / € is given, the brackets denote the duality between V and 
V . More precisely, we consider the limit boundary partial differential 
system associated with l|5.ip when e — Q. This is the membrane system, 
which according to proposition 14.61 is eUiptic, satisfies the SL on Fq 
but does not on Fi. 

5.2 Sensitive character 

Let us now recall the definition of sensitive problem. For a more 
complete description, see [7] and [15]. Let us comment a little on pro- 
position 14.61 

The SL condition is not satisfied on a free boundary when e = 
for the variational problem l|5.ip . Specifically, the membrane problem 
is of total order four for elliptic surfaces. The number of boundary 
conditions should be two. On a fixed boundary Fq they are : 

Ui M2 = 0. (5.2) 

Note that the trace of U3 does not make sense in the membane fra- 
mework. The previous boundary conditions satisfy the SL condition. 
Oppositely, on the free boundary Fi the conditions are : 

T'^f^{u)ni3 = 0. (5.3) 
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Let us admit that (|4.14p has (in some sense) a solution u. Replacing it 
in the three equations l|4.14p and in the boundary conditions on Fi of 
l|4.14p . one obtains that the corresponding T°'^{u) satisfy the elliptic 
membrane tensions system with Cauchy conditions on the part of Fi 
of the boundary. As this last problem has in general no solution in 
usual spaces, it follows that the membrane problem (|4.14p cannot (in 
general) have solution in usual spaces. We shall see that existence of 
the solution (as well as the convergence for e ^ 0) only holds in very 
abstract spaces (out of the distribution space). 

On the other hand, the boundary condition l|5.2p constitutes the 
Cauchy condition for the rigidity system "fap{u) — 0. According to the 
uniqueness theorem for elliptic Cauchy problem ( see proof of Lemma 
I4.2p an elliptic shell is inhibited (or geometrically rigid) provided that 
it is fixed (or clamped) on a part (or the whole) of the boundary. When 
the boundary is everywhere free, the shell is not inhibited. Coming back 
to the inhibited elliptic shells, we see that when the whole boundary 
is fixed, the membrane problem is classical (the boundary condition 
satisfies the SL condition). But, when a part of the boundary Fq is fixed 
whereas another one Fi is not, the boundary conditions satisfy the SL 
condition on Fq but not on Fi. This problem is out of the classical 
theory of eUiptic boundary value problems and is called sensitive for 
reason which will be self evident later. 

Let us consider formally the variationnal formulation of the mem- 
brane problem l|4.14p (i.e. with e — 0) : 



where Va is the completion of the "Koiter space" V with the norm 



The fact that [It'll a is a norm on V follows from lemma [421 
At the present state, it should be noticed that the previous comple- 
tion process is somewhat abstract and the elements of Va are not neces- 
sarly distributions. Indeed, as the SL condition is not satisfied on Fi, 
we may construct corresponding solutions with u ^ Q and "fap{u) — 
which are rapidly oscillating along Fi and exponentially decreasing in- 
wards Q.. This is only concerned with the higher order terms. When 
taking into account lower order terms (which are "small" for rapidly 
oscillating solutions), we see that we may have "large u" with "small 
7a/3(w)" (i.e. small S„,/3||7a/3(u)||L2) and then small membrane energy. 
Accordingly, the dual where / must be taken for l|5.4p to make sense 
is "very small". 

The above property originates the term "sensitive". The problem is 
unstable. Very small and smooth variations of / (even in 'D{Q)) induce 
modifications of the solution which are large and singular (out of the 
distribution space). 

5.3 Abstract convergence results as e ^ 

In this section we recall some abstract convergence results (in the 
norm of the specified spaces), see [3j and f7| for more details. 



Find u £Va such that, Vw G Va 
a{u,v) = {f,v), 



(5.4) 
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Recalling the problem we are studying, we know that the shell is 
geometrically rigid : 

V eV and -fapiv) = =^ v = 0. (5.5) 

Let A and B be the continuous operators from V into V' associated 
with the forms a and b by : 

{Au,v) — a{u,v) and {Bu,v) = b{u,v) Vu,w S V, (5.6) 

so that equation (|5.ip becomes : 

Aw^ + e^Bu^ = /. (5.7) 

Lemma 5.1. The operator A is injective and its range, TZ{A), is dense 
in v. 

The proof is not difficult, see |7j if necessary. 

It then appears that the operator ^ is a one-to-one mapping of V 
onto TZ{A), which is a dense subset of V . Let us define a new norm by 

\\v\\v,^\\Av\\v'. (5.8) 

Obviously V is not complete for the previous norm. But A defines 
an isomorphism between V (with the norm Va) and TZiA) (with the 
norm V). Automatically, A has an extension by continuity which is an 
isomorphism between the completions of both spaces. Denoting by A 
the extended operator and by Va the completion of V with the norm 
l|5.8p , A is an isomorphism between Va and V (which is the completion 
of Ti-iA) with the norm of V). Equation (|5.7p may be written as well : 

Au^ + e^Bu" = /. (5.9) 

Remark 26. In order to pass to the limit as e ^ 0, the classical way 
consists in obtaining an a priori energy estimate of by taking the 
duality product of l|5.9p with u^. But such a way needs a hypothesis 
of boundedness of the functional / with respect to the limit form a 
and this does not work for any / € V . In the general case, following 
an idea of Caillerie |3], see also [7], which consists in proving that the 
term e'^Bu^ tends to zero in V , one can pass this latter term to the 
right-hand side, and show that it tends to / in T^'. Then using the fact 
that A is an isomorphism, it is possible to prove the existence of a limit 
of in Va- Specifically we have the following result. 

Theorem 5.2. There exists a unique element u° in Va such that 

Au° = f. (5.10) 
Moreover the following strong convergence holds in Va : 

as 0, (5.11) 
where u"^ ^ V is the solution of f5.9]} . 
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The proof, which follows the trends outlined above, may be seen in 

m 

Remark 27. It should be emphasised that theorem 15.21 holds true wi- 
thout special hypothesis on / (besides the obvious one / G V). The 
Hmit M° S Va is the solution of the abstract problem (|5.4p , which is not 
a variational one. The classical variational theory of the Hmit needs a 
supplementary hypothesis on / : there exists C > such that 

\\{f,v)\\<Ca{v,vy/\ VveV, (5.12) 

which is very restrictive in shell theory. 

For the sake of completness, let us give the elements of the classical 
Hmit theory under the assumption (|5.12p . 

We first note that in such a case, a{v,vy^'^ defines a norm on V. 
Let Va be the completion of V with respect to that norm (which should 
not be confused with Va)- We then note that l|5.12p shows that / may 
be extended by continuity to an element of V„'. We shall denote this 
extension by / again. Obviously, the variational problem 

Find M° e Va such that, Vu € K ik -lo^ 

a(uO,«) = (/,«), 

is well posed and has a unique solution. We then have the classical 
convergence result (see [10] e.g. or even |18j l 
Theorem 5.3. Under the assumption f5.1S^} . we have 

M° strongly inVa as e ^ 0, (5-14) 

where and m° are the solutions of ^5.1]) and 115. 10\) respectively. 

Let us now briefly recaH the non-inhibited case when l|5.5p does not 
hold. In such a situation, there is a convergence result towards a limit 
with vanishing membrane energy. More precisely, we define the kernel 
G of a : 

G^{veV; -/afiiv) ^0}^ {v e V; a(w, v) = 0}. (5.15) 

It is to be noticed that G is a Hilbert space with the norm of V. But 
as a{v,v) = in G, we see that the norm of y in G is equivalent to 
5(w,i;)^/^. As a consequence, the problem 

Find -y" e G such that, Vw S G -ifW 

biv',w) = {f,w), ^^■^^> 

is well posed and has a unique solution. Moreover, since the "Hmit 
form" a in (|5.ip vanishes on G, it implies some kind of weakness in G. 
The solution will be very large and we should define a new scaling in 
order to have a finite limit, w"^ = e^u*^ , l|5.ip becomes 

Find G V such that, \fw G V ik i7\ 

e^'^a{v'^ ,w) + b{v'^ ,w) — {f,w), ^ ' ' 

we then have, see [16] e.g. for the proof 
Theorem 5.4. Under the assumption G 7^ 0, 

w° strongly in V, (5.18) 

where and w° are the solutions of ^5.17]) and f5.16]) respectively. 
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6 Heuristic asymptotics in the previous pro- 
blem 



The aim of this section is the construction, in a heuristic way, of an 
approximate description of the solutions u'^ of the Hnear Koiter model 
for small values of e. Indeed, coming back to the Koiter problem for 
£ > 0, in the sensitive case, the problem is not really to describe the 
limit problem (which in general has no solution in the distribution 
space; in particular the space Va (see l|5.8p ) where there is always a 
limit, is not a distribution space), but rather to give a good description 
of the solution for very small values of e. This we shall try to do. We 
shall see that heuristic considerations allow to construct a simplified 
model accounting for the main features of the problem. 

To do so we shall use the heuristic procedure of [7]. In this latter 
article, we addressed a model problem including a variational struc- 
ture, somewhat analogous to the problem studied here, but simpler, 
as concerning an equation instead of a system. It is shown that the 
limit problem contains in particular an elliptic Cauchy problem. This 
problem was handled in both a rigorous (very abstract) framework and 
using a heuristic procedure for exhibiting the structure of the solutions 
with very small e. The main difference is that in the present work, we 
deal with systems instead of single equations. 

We shall see that heuristic considerations involving minimization 
of energy allow us to reduce the problem to another on the boundary 
Fi. In that context, it is seen that the "pathological" operator A is 
represented by a smoothing operator S (i.e. sending any distribution 
to a C°° function), whereas the "classical" operator B is represented 
by a "classical" elliptic operator Q. Denoting by s{x,^) and q{x,£) the 
corresponding symbols (here x is the arc on Fi), s is likely exponentially 
decreasing for ^ ^ oo, whereas q is algebraically growing. The action 
o{ S + e^Q on test functions is given by : 



{S^e^Q)e{x) = {2^)-^ e'«-[s(x,e)+e'g(x,C)]0~(e) de (6.1) 



It is then apparent that, when e is small, operator S is significant 
only for bounded values of ^, whereas e^Q describes the behavior for 
^ ^ oo. If 1^1 << log{l/e), then the symbol of the operator S + e^Q is 
equal to (l+o(l))s(a;, ^) and for |^| >> log{l/e), it is {l+o{l))e^q{x,S,). 
The balance of S and e^Q is obtained for values of ^ such that : 



This is the window of frequencies allowing a good description of 
the simultaneous infiuence of S and e^Q, which is precisely our aim. 
Moreover, it is easily seen that the range of frequencies l|6.2p is respon- 
sible for most of the contribution to the integral l|6.ip . This property 
is of great interest for the construction of the heuristic approximation. 
More precisely, the heuristics incorporate approximations for large |^|. 
This amounts to saying that only the most singular parts of the so- 
lutions are retained, or equivalently, that the approximate solutions 




\^\^log{l/e). 



(6.2) 
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are defined up to more regular terms. This is for instance the kind of 
approximation which is used in the construction of a parametrix. We 
also note that, as l|6.2p involves "moderately large" values of |^|, the 
"general quality" of the approximation is not very good, as it is only 
accurate for very very small values of e. 

It should be noticed that numerical computations [2j carried out 
with very rehable software (including an adapted mesh procedure) for 
the Koiter problem with very small values of e agree with the ove- 
rall trends of our heuristic procedure. It appears that most of the de- 
formation consists in very large deformations along Fi exponentially 
decreasing inwards (then in good agreement with the "local lack 
of uniqueness" implied by the non-satisfied SL condition). As e de- 
creases, the amplitude increases, whereas the wave length decreases 
very slowly, verifying fairly well (|6.2p . The paper p] also contains nu- 
merical comparisons with the case when the shell is fixed all along its 
boundary, which is classical (as the SL condition is satisfied all along 
the boundary). The differences are drastic for small values of e. 

6.1 Introduction to the heuristic asymptotic 

A first remark in the context described above is that sensitive pro- 
blems may be considered as "intermediate" between "inhibited" and 
"non- inhibited". Indeed, "inhibited" means that v G V and 70/3(1;) = 
implies v = 0, whereas "non-inhibited" means that there are non vani- 
shing elements v ofV such that 7a/3(w) = 0. Strictly speaking, sensitive 
problems enter in the class "inhibited", but there are non vanishing 
elements vofV with "very small" 70,^3 (w). 

In order to minimize the energy 

a{v,v) + e^biv,v)-2{f,v), (6.3) 

it is clear that we may proceed as in non-inhibited problems. The 
solution with small e "avoids" the (larger) membrane energy a, so that 
roughly speaking, solutions for small e should have 7q/3(w) vanishing 
or at least very small with respect to v. 

Obviously, it is impossible to impose the four boundary conditions 
l|4.23p on Tq with the "exact" system jafiiv) = as they imply v = 0. 

Nevertheless, we shall see in Section 16.2.11 that it is possible to 
construct functions satisfying the two boundary conditions Un — ut — 
on Tq with the "non exact" system 7q/3(w) = in the sense that japiv) 
will be "very small" (i.e. SQ,^^||7Q/3(w)||i2 will be very small). This will 
imply a "membrane boundary layer" in the vicinity of Tq involving the 
bilinear form a. To this end, we shall first construct a set of functions 
V with only one vanishing component on Tq. Choosing (for instance) 
the normal component, we define : 

G° = {v, jaisiv) = on 17, i;2 = on Tq}, (6.4) 

the regularity is not precised as we shall later take the completion, we 
may consider C°° functions for instance. It is to be noticed that w is a 
triplet of functions. 
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Recalling Remark[22l we know that up to a finite dimensional space 
composed of smooth functions, the space G" is isomorphic to the space 
of traces on Fi : 



In the sequel, when we will consider a function w G G", we will 
consider a function of the equivalence class for the quotient operation 
described in Remark [22l Moreover, we shall consider indifferently the 
functions w obtained after a quotient operation on fl (for the finite 
dimensional space) or their traces w on Fi. 

Moreover, the conditions U3 — dnU^ = on Fq of (|4.23p will be 
satisfied with the help of a "fiection sublayer" involving the bilinear 
form b ; its effect is not relevant (see Section I6.2.2p . 

According to the previous considerations, we shall consider the mi- 
nimization problem on G° instead of on V . This modified problem ob- 
viously involves the a-energy and the e^6-energy. A natural space for 
handling it should be the completion G of G" with the corresponding 
norm. 

The fact that we may "neglect" the functions in the finite dimension 
space of smooth functions follows from the fact that we are interested 
in the singular part. 

6.2 The boundary layer on Fq 

Let w be in G'^ (see (|6.4p ) and let e > be fixed. The aim of this 
section is to build a modified function zw" of w in a narrow boundary 
layer of Fq in order to satisfy the supplementary boundary conditions 

= = dnW^ = on Fq. 

The present problem is analogous to the "model problem" of [7] in 
the case of a singular perturbation, i.e. U Section 7.1.2. Indeed, the 
membrane problem is of total order 4 allowing 2 boundary conditions 
{wt — Wn — on Fq, whereas the complete Koiter shell problem is of 
order 8, allowing 4 boundary conditions (we shall add W3 = dnW^ = 0) 
on Fq. It appears that the two first conditions {wt = w„ = 0) may 
be obtained from elements of by modifying them on account of a 
"membrane layer" which relies on the membrane system, of thickness 
of order iog(\yg-) on Fq, whereas an irrelevant boundary layer will be 
considered in Section fG. 2. 21 

6.2.1 The membrane boundary layer on Fq 

In this subsection, we proceed to modify the element w of G° in 
order to satisfy both conditions wi = U2 = on Fq. 

Let Fq be a neighborhood of Fq in disjoint with Fi and suffi- 
ciently narrow to be described by the curvilinear coordinates yi — arc 




(6.5) 




(6.6) 
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of To and y2 = distance along the normal to Fq. Let {■i/jj{yi))j^j be a 
partition of the unity associated with Fq and let rj e C°°(M+;M+) be 
a cut-off function equal to 1 for small values of y2. 

The mappings 9j defined by 9j{yi,y2) = ^]{yi)v{V2), where y2 is 
the (inwards) normal coordinate along Fq , define a partition of unity 
in Fq ; in particular, for a given G G"^, we have : 

V(2/i,?;2) e fo, w{yi,y2) = T:jejOj{yi,y2)w{yi,y2)- (6.7) 

Let us now fix j in J and y2 such that (2/1,2/2) G Fq, the function 
^j('i y2)w'(', 2/2) has a compact support, we denote by w^{-,y2) its ex- 
tension by zero to R and by J-{w^) the tangential Fourier transform, 
2/1 ^ Ci , of w-J . 

Let us first exhibit the local structure of the Fourier transform of 
close to Fq . Denoting by 9j the multiplication operator by 9j , recaUing 
that the commutator of the operator 7 associated with 7^/3 and 9j, 
denoted by [7, 9j\, is a differential operator of lower order, taking the 
7 operator in the new coordinates (2/1,2/2) (which, according to our 
approximation close to Fq, has the same principal part) and using that 
w € G", we see that : 



for some i > 0, Uap being differential operator of order less than the 
order of 7^/3 . 

Now, according to the general trends of our boundary layer approxi- 
mation, we can neglect the terms of lower order in l|6.8p and we can 
proceed as in the construction of a parametrix (freezing coefficients, 
dropping lower order terms, solving such simpler equation via tangent 
Fourier transform and gluing together the solutions for different j), so 
that l|6.8p becomes 



and taking the tangential Fourier transform denoted by JF(w-')(^i, 2/2) 
this yields 



lc.p{w') + Uc.fi[y,D)w' = on K X (0,i), 



(6.8) 



7„^(wJ) = on M X (0,<). 



(6.9) 



The previous system can be rewritten as 




(6.10) 




(6.11) 



with 




The general solution of the system l|6.1ip is : 
JP(w^)(Ci,2/2) = v4w+e^+(«i)^^ +Bw-e^ 



,A-(«i)y2 



(6.12) 
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1 |,A+(e 

\ 622 

tVbiib22-bl, and A_(a) = "^Cife " ^T^il^ 

Since w G Gq, it follows that T{w2){£.i,0) = hence A 




"12 • 



-B. 



Consequently, we deduce that 
&11622 



^2 
"12 



(^i)(6,o)( 



(6.13) 

This expression exhibits the structure of (the Fourier transform of) 
w in a narrow neigbourhood of Tq. It was expressed in terms of (the 
Fourier transform of) the trace of its third component on Tq, but this 
choice is arbitrary. 

We now proceed to the modification of in w^"" in a narrow boun- 
dary layer of Fq in order to satisfy (always within our approximation) 
the equation coming from l|4.22p for £ = (this is the membrane boun- 
dary layer associated with the membrane system of Section [473| . Using 
considerations similar to those leading to l|6.8p . this amounts to 



(6.14) 



(7M17) w^'' + U{y, D)w''' = on M X (0, t), 



where J7 is a differential operator of lower order than four, 7* denotes 
the operator : 



7 



and 



A 



A 





- 


bii \ 
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b22 


d2 


di - 


26i2 / 


1111 


^1122 


^1112 


2211 


^2222 


^2212 


1211 


^1222 


^1212 



Therefore dropping as before terms of lower order, we have : 



(7*ii7)?I)^'' ==OonRx (0,t), 
which can be rewritten as 

f(7o -71*32)^1(70 +71^2) l^i-'" = on M X (0,i), 



— T 

with 7* = 7 and 



(6.15) 



(6.16) 
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Hence taking the tangential Fourier transform, we look for solutions of 
the system : 



d 

d?/2 
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with 









At this moment, it is worthwhile to compare l|6.17p and l|6.1ip . We 
see that the given function (rather its Fourier transform) solves the 
"right half" of (|6.17p . i.e. the expression on the right of Ai in (|6.17p . 
Obviously, the "left half" accounts for the "adjoint part", coming with 
integration by parts from the bilinear form a (see (|2.5p ). Our aim in 
constructing the modified w^"" is to satisfy the conditions = iDj'^ = 
on 2/2 = 0, whereas for "large j/2" (in the sense of "out of the layer") 
the modified w-'" coincides (up to small terms) with the given . We 
now proceed to write down the general solution of (|6.17p on account 
of its special structure. 

For A e {A_(^i), A+(^2)}, let us consider the function k defined by : 

kiiuy2)^{y2W + v)e^y', (6.18) 
where w e {w_,w+} is a solution of 

(jo + X^i^w ^ 0, 

and V is unknown. We then search for solutions of l|6.17p under the 
form l|6.18p i.e. : 

{(^'^ ~'fi^^)M% + ii^j)Hii,y2) = o, (6.19) 



(7o+7i-^)(2/2W + e^^'w) = ((7o + A7i)?; + 7iw)e^^^ 



We check that 

d 

So that (|6.19p becomes 

((70^ - 7^ "d^)^i('^o + '^^"d^O (^^^ ^ 

dy2 ^ ^ 
This amounts to saying that Ai ^(7o + A7i)w+7iw^ is an eigenvector of 

2"' / T' 

7o — associated with the zero eigenvalue. Since dim Ker (70 — 

X-Ji^ — 1, denoting by uq a non vanishing vector of Ker ^70^ — At-^ ^ , 
then V should satisfy 

(70 + A7i)w + 7iw — A^^{tuo), for some t e C. (6.20) 

According to the Fredholm alternative, a necessary and sufficient condi- 
tion for the existence of such a f is that 

A^^{tuo) — 71 w e (Vect uq)'^. 
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Since Ai is positive definite, we deduce that {A^ uq, uq) > hence 
T — (71". "o) satisfies 

(Aj^ «o,Mo) 

{tA^^uo - ^iw,uo) = 0. 

It follows that the vector v £ exists and is unique (up to an 
additive and arbitrary eigenvector, which is irrelevant in the sequel). 
Consequently, k defined as above satisfies (|6.19p . 

Repeating this argument twice (first for A+(^i), and then for A„(^i)), 
and denoting by w+ and w_ the corresponding vectors v, we see that 



with arbitrary Ci, C2, C3, C4 is the general solution of i|6.17p . 

We are now determining Ci, C2, C3, C4 in order to satisfy the boun- 
dary conditions = d2wl°' = at ?;2 = and the "matching condi- 
tion" with , i.e. in the context of boundary layer theory (for large 
l^il), w^" should become out of the layer. 

Let us now explain the process of matching the layer : out of the 
layer, we want w^'^ to match with the given function Since |^i| >> 

1, then |Ci|y2 >> 1 and ^''^^j^^^ —\S,i\y2 >> 1 which means that 
1/2 » I o TTT (but we still impose that 2/2 is small in order 

to be in a narrow layer of Fq where (|6.13p holds) ; this is perfectly 
consistent, as we will only use the functions for large |^i|, hence the 
terms with coefficients C2 and C4 are "boundary layer terms" going to 

zero out of the layer (i.e. for |?/2| >> ^(]|7j))- 

The matching with l|6.13p out of the layer then gives 



C3 = and Ci = 
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•^(*'3)(a,o). 



2|eilV^ll&22-&?2 

The two other constants C2 and C4 are determined by 

.F(ii^")i(a,0) = and T{w^'')2{iu^) - 0, 
which yields the existence of two constants a and /3 such that 

C2 = aCi and C4 = /3Ci. 
So that the modified solution is of the form : 



3.22) 



'22 



>^+{(,i)y2 



2|6lv/^ii622-6?2 
+ ((a + (3v2)w- + /?^;_)e^-(«i)^^).F(w|)(a,0) 



.23) 



The modification of the function Wj then consists in adding to it 
the inverse Fourier transform of 

^'11622 



"12 



((a + 1 + /3y2)tl!- + /3w_)e^-«i)^^.F(zZ;^)(6, 0). 

(6.24) 



,^+(?l)y2 



(6.21) 
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We shall study in the sequel the behavior of such an expression. The 
role of the constants a and /? is not relevant, and we may assume, for 
instance that a = — 1 and /3 = 1 (this amounts to change w_ and -5-). 
As the result, the modification of the function Wj consists in adding to 
it the inverse Fourier transform of 



22 



= (2/2^_ +t.^)e^-(«^)^^.F(u;^)(ei,0). (6.25) 



2|eilV^ii622-fo?: 

More precisely, on account of considerations at the beginning of 
Section 6 (see in particular l|6.ip and (|6.2p ). the modification should 
only be effective for large |^i|, accounting for "singular parts" of the 
solution. Moreover, in order to have w°' G V ^ we shall also impose 
= d2W-'^ = on To (the other two conditions = 9„wf = 
on To will be adressed in Section [6. 2. 2p . To this end, we multiply the 
added term by a cutoff function avoiding low frequencies (It should be 
remembered that this is one of the typical devices in the construction 
of a parametrix) . More precisely, on account of l|6.2p , we shall only keep 
frequencies of order more or equal than [Zog(l/e)]^/^, which preserve 
the useful region (|6.2p and are large (then consistent with the fact that 
the modification is a layer). Moreover, in order to the modified function 
satisfy the boundary conditions, we must also take into account the low 
frequencies of the Fourier transform which we multiply by a smooth 
vector p{y2) such that pi(0) = P2(0) = and p{y2) = for yi > C 
for a certain C. The division into high and low frequencies is defined 
by a smooth function H{z) equal to 1 for > 1 and vanishing for 
l^l < 1/2, with z = [;og(i/g)]i/2 ■ Finally, we define the function 

Me,e,y.)^(l-i^(j^^;^(0^ (6.26) 

216 1 Vbnh2 b'i, V^'^- ^"-)' ^ [logil/eWf^ ^' 

which obviously has its first and second components vanishing for j/2 = 
0. Now we can modify the function Wj by 

Swj =Wj -Wj, (6.27) 

where dwj is defined by its Fourier transform : 

t(^Sw,) ^Tiwi)i^i,0)h{e,^,y2). (6.28) 

Remark 28. The constant C in the definition of piyi) is chosen suffi- 
ciently small for this function to vanish out of the layer of close to 
Fq where the curvilinear coordinates yi,y2 operate. Rigorously spea- 
king, the rest of the expression should also be multipHed by a cut-off 
function vanishing for y2 > C, but this is practically not necessary, as 
this part is exponentially small for large 

Hence summing over j and defining on Fq the family (with para- 
meter ?/2) of pseudo-differential smoothing operators (5cr(e, Di , j/2) with 
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symbol : 

<Sa(e,a,y2) = , '^'^ ,2 (y2W-+v\e^-(^^^y\ (6.29) 

2|Cl|V^ll&22 - 6l2 ^ ^ 

we see that the modification of the function w : 

Sw^w" -w (6.30) 

is precisely the action of Sa{e, ^2) on 11)3(2/1, 0). 

Once w"" is constructed, it is worthwhile computing its a-energy. 
This we proceed to do. More generally, we shall compute the form a 
for two functions u° and w°-. 

Let us now compute the leading terms of the a-energy of the mo- 
dified function w'^ . 

Let V and w be two elements in and v",?!)" the corresponding 
elements modified in the boundary layer. As the given v and w satisfy 
la/siv) = Ja/iiw) = 0, the a-form is only concerned with the modifica- 
tion terms 6v and 6w. Then, within our approximation, we have : 

+00 



a(«^w")=/ A'^f'^^dyi JM^ihx^.iSw) dy2. (6.31) 



where the integral in dy2 is only effective in the narrow layer. Using 
the partition of the unity 9j and denoting as before by Swj{-,y2) the 
extension with value to E of dj{-,y2)Sw{-,y2), we have 



a(^)^u.'^) = E,■fc / A'^'^'^dyi lo.^i6v,)-fx^,iSwk) dy2. (6.32) 



Consequently, using the tangential Fourier transform yi £,1 and 
the Parceval-Plancherel theorem, dropping lower order terms (within 
our approximation, we only consider expressions with large |^i| which 
amounts to take H = I in (|6.26p ). we deduce that 

a(w°,i&'') = 

E,, y ^i^^i/, (7o+7i^)Me,e,2/2)^(«^)(Ci,0)) X 



70 + 71 ^) S<y{e, t y2)^{wlmi , 0)) dys 



A, d^i / , " :((7o + A_7i)t'-+7i^-)e^-^-.F(i)i)(6,0)) x 

-00 Jq 2|^i|\/6ii&22 - fcl2 ^ ^ 



bub: 
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2|Cllv/^ll&22-fo?2 



Hence, on account of the definitions of 70, 71, A_ and W- integrating 
in 2/2, we know that 



a(«",i&'') = S.-fc / 0\£i\ni')3\v.=onw'')3\v,=oh'{e,£,y2) d^i, 

(6.33) 
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with e = 6'(A"^^^, (w_)i(0),6a/3,/i_), where /i- = is independent 
of^i- 

This expression (|6.33p only depends on the trace {v^)^y^=Q{yi) and 
(w'^)3|y2=o(yi)) which are functions defined on Fq. 
Remark 29. The important fact in (|6.33p is the presence of |^i|. This 
comes from /q^°° e^^-^^ dj/2 and analogous, on account that this inte- 
gral is equal to 

We now simplify this last expression using a sesquilinear form in- 
volving pseudo-differential operators. 

Then, defining the elliptic pseudo-differential operator P{yi,-^) 
of order 1/2 with principal symbol 

{0\^i\y/^h{e,ty2), (6.34) 
and summing over j and k, we obtain 

a{v\w'^) = ^ P(^)(^)3)|^^p(^)(^3)|^^^ ds. (6.35) 

Remark 30. As we only considered the principal terms for large |^i|, 
we may define as well P{£,i) by the symbol 

P(6) = ^(l + ICiP)'/'- (6.36) 

The corresponding pseudo-differential operator is elliptic of order 1/2. 
Remark 31. We shall use the definition l|6.36p . which is more pleasant 
than (|6.34p . as such a P defines an isomorphism from H''(ro) onto 

i7'*+i/2(ro), s e M. 



6.2.2 The flection sublayer on Fq 

The structure of the flection sublayer, see the beginning of Section 
6, accounting for the two new boundary conditions — dnW3 — 
follows from classical issues in singular perturbation theory, as in [7] 
section 7.1.2, [22j and It is mainly concerned with a drastic change 
of the normal component ws (whereas the conditions on wi and W2 are 
satisfled). The speciflc structure is analogous to the layer in [19]. 

The thickness is of order S = e^l"^ . This may be easily seen by 
taking into account only higher order terms in the membrane and the 
flection systems ; eliminating w\ and W2 , we obtain an equation for wz ■ 
The membrane terms are of order 4 and the flection terms are of order 
8. In the layer, the derivatives of order n have an order of magnitude 
0(1^). As both membrane and flection terms are of the same order of 
magnitude in the layer, we thus have 

O(^) =e O(^), 

which furnishes 6 — ©(e^/^). 

It is easily seen (as in [7] section 7.1.2) that the presence of this 
flection sublayer plays a negligible role in the asymptotic behavior. 
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Indeed, proceeding as in the previous membrane layer, we see that the 
expression analogous to (|6.35p has the form : 

e^aoiv^w^) = e' Po(^)(^5)|ro^o(^)(^)|ro ds, (6.37) 

where Pq is an operator of order 0. Going on to next Section 16. 4^ 
the action of sublayer amounts to change ^ to ^ + e^C where C is a 
smoothing operator. Equivalently, we may change B to B + C which is 
again a 3-order operator (as C is smoothing) . The asymptotic behavior 
does not change. Equivalently, in l|6.ip . the effect of the sublayer is to 
change s to s + e^so where sq is a smoothing symbol, or g to g + sq 
which is again the symbol of an operator of order 2m > 0. 

For that reasons, the influence of the sublayer will no more be 
mentioned. 

6.3 Formulation of the problem in the heuristic 
asymptotics 

Presently, our aim is to formulate problem l|5.ip on the space of the 
u° with u G G°. The forms b{u,v) and (/, f) should be written in the 
framework of our formal asymptotics, for and obtained from u 
and V defined on Fi by solving (|6.6p and modifying u and v with the 
Fo -layer. 

The computation of the 6-energy form is exactly analogous to that 
of [7] Sec. 5.3. It follows the ideas of the previous section in a much 
simpler situation. As only the third component is involved in the higher 
order terms of the form b (see l|2.15p and (|2.2p ). we have 

6(u^^)'^) « / B'^l'^^^dapu^^dx^i^ d^ dx. (6.38) 

Moreover, from (|6.4p - (|6.6p and according to our approximations 
analogous to the construction of a parametrix, u, v are only signifi- 
cant in a narrow layer adjacent to Fi. The local structure is analogous 
to (|6.12p where obviously the decreasing solution inwards the domain 
should be chosen. This gives the obvious local asymptotics 

^3(e,^/)=«3(a)e^-'«^'^^ (6.39) 

where A_(fi) is proportional to |^i|. After substitution (|6.39p in l|6.38p 
a computation analogous to that of Section 16.2.11 (but much easier) 
gives (using a partition of unity) : 

/ + 00 
Ofe (yi) 161^3(6)^3^ (6 )d6 
-oo 

where (jkivi) are smooth positive functions on Fi depending on the 
coefficients. The function |6|^ comes obviously from the integrals in the 
normal direction of products of second order derivatives of functions 
of the form e'^-'^^^^^, with A_(6) proportional to |<^i|. 



34 



Then, defining the pseudo-differential operator Q(gf-) of order 3/2 



with principal symbol 



v/C(yi)iaP, 

we have within our approximation : 



(6.40) 



f B"^^'^do^pUsdxf,V3dx= f Q(J-)u Q{^)v dyi. (6.41) 
Jn Jr, dyi dyi 

We observe that the operator Q is only concerned with the trace 
on Fi and yi which denotes its arc. 

The formal asymptotic problem becomes : 



Find e G such that Vw G G 



where G is the completion of G" for the norm 



(6.42) 



ro 



on 



Q{v:i) 



ds 



Remark 32. For e > 0, (|6.42p is a classical Lax-Milgram problem. 
Continuity and coerciveness follow from the ellipticity of the operators 
P and Q. 



6.4 The formal asymptotics and its sensitive beha- 
viour 

In the sequel, we shall denote 

a{v',w) = ^ Pi^)P(^)ds (6.43) 

/3('D^w) = / Q[v') Q{w) ds. (6.44) 

We observe that the problem l|6.42p is again in the same abstract 
framework as the initial problem l|2.17p . Nevertheless, the context is 
different, as the non-local character of the new problem is apparent 
from the structure of the space G. Let us define the operators 

Ae£{G,G'), BeC{G,G') (6.45) 

by 

a{v,w) — {Av,w) (3{v,w) — {Bv,w). (6.46) 
Let G^ be the completion of G with the norm 

\\v\U = \\Av\\g'. (6.47) 

Denoting again by A its extension to £(G^,G'), which is an iso- 
morphism, we may rewrite (|6.42p in the form : 

[A + eB^i' ^ F, (6.48) 
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where F e G" is defined by 

{F, w) ^ fw dx, Vw e V. (6.49) 
Jn 

It follows that 

-C^ — > 5" strongly in Ga, (6.50) 

where 

= F. (6.51) 

Reduction to a problem on Ti 

In order to exhibit more clearly the unusual character of the pro- 
blem, we shall now write l|6.42p in another, equivalent form involving 
only the traces on Fi . Coming back to l|6.6p . let us define TZq as follows. 
For a given w G C°°(ri) we solve (|6.6p and we obtain 

= TZqw. (6.52) 

Using the regularity properties of the solution of (|6.6p . it follows 
that TZi)W is in C°°(ro). Moreover, we may take in (|6.6p a w in any 
H'^iTi), s e M and the corresponding solution is of class C°° on Fq 
and its neighbourhood, so that TZq has an extension which is continuous 
from WiTi) to C°°(ro). We shall denote by TZq such an extension, so 
that 

7^o e C{H'{Ti),H-{To))ys,r £ M. (6.53) 
Then, l|6.42p may be written as a problem for the traces on Fi : 

f Find G H^' ^{Ti) such that Vu; G H^/^{Ti) 

i /ro P(i)7^o^;^P(i)7^o^« ds + £2 /^^ g(^)„^ Q(i)u; ds ^ dx, 

(6.54) 

where the configuration space is obviously H^/'^{Ti). The left hand side 
with e > is continuous and coercive. 

Remark 33. Coerciveness follows from the ellipticity of Q; as it is of 
order 3/2. Strictly speaking, this only ensures coerciveness on the lea- 
ding order terms, which may "forget" a finite-dimensional kernel. But 
this is controlled by TZq, as it is a surjective operator. Indeed, TZ^v ~ 
impHes 7a/3(w) = with = V2 = Q on Tq, which implies v — Q 
(and then d = 0) using the uniqueness of the Cauchy problem for the 
rigidity system. 

Here F G H-^/^{Ti) is defined for / G by 

(-P',w>i/-3/2(ri),_f/3/2(ri) = (6.55) 

We note that, for instance, when the "loading" / is defined by a 
"force" F on Fi, this function is the F in (|6.54p . Obviously, (|6.54p may 
be written : 

From l|6.53p we see that T^q ^l^o a smoothing operator, i. e. : 

ni G /:(i?-'^(Fi),i7-^(Fo)),Vs,r G R. (6.57) 
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Now we define the new operators (but we use the same notations) 

A = n*„P*Pna e C{H'{Ti),H''{To)),y3,r eR, (6.58) 

B = Q*Q e/:(i/3/2(ri),i7-3/2(r,)). (6.59) 

Obviously, B is an elliptic pseudo-differential operator of order 3, whe- 
reas ^ is a smoothing (non local) operator. Then l|6.56p becomes 

(A + e^By" = F in H-^^^iVi). (6.60) 

Once more, the problem (|6.54p is in the general framework of (|2.17p . 
so that we can define the space V = H^/^(ri) and its completion Va 
with the norm 

llwlU = ||^w||H-3/2(ri)- (6.61) 

Denoting similarly by A the continuous extension of A, which is an 
isomorphism between Va and V, we obtain 

u" strongly in Va, (6.62) 

where £ Va satisfies 

Au'^ = F. (6.63) 

Obviously, this equation is uniquely solvable in for F £ V — 
H~^/^(ri). But, the unusual character of this equation appears now 
clearly : 

Proposition 6.1. Let F G iJ-3/2(ri) and F ^ C°°(ri), then the 
problem ^6. 63]) has no vP solution in I?'(ri). 

Demonstration. If vP £ V'lT^^) was a solution of (|6.63p . as Fi is com- 
pact, vP should be in some T?'*, then recalling (|6.58p . we should have 
AvP e C°°(ro), which is not possible. Moreover, (|6.60p is clearly of the 
form l(6TT|) . □ 
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